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ABSTRACT 



We present spectroscopic measurements of the [0 ill] emission line from two 
subregions of strong Lya emissi on in a radio-quiet Lyman-alpha blob (LAB). 
The blob under study is LABI jsteidel et all boooh at z ~ 3.1, and the [Om] 
detections are from the two Lyman break galaxies embedded in the blob halo. 
The [O ill] measurements were made with LUCIFER on the 8.4m Large Binocular 
Telescope and NIRSPEC on 10m Keck Telescope. Com paring the redshift o f 
the [O ill] measurements to Lyct redshifts from SAURON (jWeijmans et al.ll2010i ) 
allows us to take a step towards understanding the kinematics of the gas in the 
blob. Using both LUCIFER and NIRSPEC we find velocity offsets between the 
[O ill] and Lya redshifts that are modestly negative or consistent with km 
s _1 in both subregions studied (ranging from -72 ± 42 - +6 ± 33 km s _1 ). 
A negative offset means Lya is blueshifted with respect to [O ill] a positive 
offset then implies Lya is redshifted with respect to [Om]. These results may 
imply that outflows are not primarily responsible for Lyman alpha escape in this 
LAB, since outflows ar e generally expected to produce a positive velocity offset 
(IMcLinden et al.ll201ll ). In addition, we present an [Om] line flux upper limit 
on a third region of LABI, a region that is unassociated with any underlying 
galaxy. We find that the [O ill] upper limit from the galaxy-unassociated region 
of the blob is at least 1.4 - 2.5 times fainter than the [O ill] flux from one of the 
LBG-associated regions and has an [Om] to Lya ratio measured at least 1.9 - 
3.4 times smaller than the same ratio measured from one of the LBGs. 
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Subject headings: galaxies: high redshift — intergalactic medium — galaxies: 
kinematics and dynamics 



INTRODUCTION 



Lyman-alpha (Lya) first became a useful tool for obser ying high-z sources with the 



discovery of large s amples of Lya emitting galaxies (LAEs) ( ICowie fc Hul Il998l ; iHu et al. 



19981 : lRhoadsll2000l ). The same narrowband imaging techniques that uncovered LAEs be- 
gan uncovering a different set of objects that were also very bright in Lya. These rarer, 
more extended, and more luminous objects are what we now call Lya blobs (LABs) (e.g., 



Steidel et al. 



2006 



2000 



Francis et al. 2001 



Prescott et al. 



Matsuda et all 12004 : iDev et all 120051 : iNilsson et al. 



20121 ) . LABs are extremely large (~ 30-200 kpc) radio-quiet Lya neb 



ulae in the high redshift universe. LABs are highly luminous (L>Lyc 



10 



43-44 



ergs s 



_1 ), and 



yet despite rigorous study in the last decade, the mechanism(s) that power this immense Lya 
flux is not fully understood. This paper will focus on investigating the kinematics of and 
me chanisms powering such objects by investigating LABI, a z ~ 3.1 LAB first discovered 
bv lsteidel etaD (boooh . 



There are currently three most widely discussed scenarios to explain both the large 
spatial extent and powerful Lya flux of these blobs. The first of thes e is that the gas i n 
LABs is heated by photoionization from massive stars and/or AGN ( jGeach et al.l 120091 ). 
A second scenario p r oposes that gas in LA Bs is heated by cooling flows / cold accretion 
( IHaiman et al.l l2000t iDijkstra fc Loebl 120091 ). Finally some authors have proposed LABs 
originate from overlapping supernova remnants from massive stars after a powerful starburst 
( jTaniguchi fc Shioyall2000l : lOhyama et al.l 120031 ) producing superwinds. 



Adding to the controversy, observations in recent years from different authors have led 
to different conclusi ons about which of these scenarios are responsible for said observations. 
Nilsson et al.l (120061 ) have argued that their observations of a z ~ 3.16 LAB were best matched 



by coo ling flows onto a dark matter halo. This is in contrast to the conclu sions of lHaves et a 



( 1201 it ), who found evidence of polarized Lya radiation in LABI. The lHayes et al.l (120111 ) 



results suggest that Lya photons are scattered at large radii from their production sites and 
this observation seems to not only favor the role of scattering in outflows in LABs, but the 
au thors contend their di scovery can actually rule out most inflow models. But a similar study 
by lPrescott et al.l ( 120111 ) found no evidence of polarization in a LAB at z ~ 2.656 and these 
authors argue their results are inconsistent with sp herical outflows an d Lya scattering from 
nearby AGN. Yet another conclusion is reached by lYang et al.l (1201 ll ) whose observations of 
Lya and Ha emission in two z ~ 2.3 LABs rule out simple infall models and models that rely 
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on large outflows. This diversity of conclusions may mean that there are diverse mechanisms 
powering different blobs (or multiple mechanisms at play in single blobs) or it may mean 
that truly conclusive observations have not yet been presented. 

To try to provide new da ta to differentiate amongst the possible LAB sources we focus 
on LABI (ISteidel et al.l 120001 ). LABI resides i n SSA22, a protoclus ter region at z ~ 3.1 



( jSteidel et al.lll998[ ). LABI e xtends ~ 100 kpc (jWeiimans et al.ll2010l ) and has a Lya lumi 



nosity of 1.1 x 10 44 erg s 1 (IMatsuda et al.l 120041) . This makes LABI one of the brightes t 



and most spatially extended LABs yet observed (ISteidel et al.l l2000l ; IMatsuda et al.l 12004 ). 
LABI is comprised of five separate regions of Lya emission known as Cll, C15, Rl, R2, and 
R3 (see Figure 1 of Weijmans et al (201 0)). Cll and C15 are both Lyman break galaxies 
(LBGs) identified b y lSteidel et al.l (120001 ). R3 has been identified as an extremely red galaxy 
(IGeach et al.l 120071 ) and may or may not be associated with a bright submillimeter source 
and nearby radio source (Chapman et al. 2001, Chapman et al. 2004, but see also Kohno 



et al. 2008, Yang et al, 20 12). Rl and R2 are not identified with galaxies (I Weijmans et al. 



20101 ). IGeach et al.l ( 120091 ) also demonstrated that AGN activity is n ot significant in LABI 
as LABI is not detected in a 400 ks Chandra exposure. We note that IMatsuda et al.l (120041 ) 
identified C15 as a separate LAB, but in this pa per we refer to all fiv e su bregions as com 



ponents of LABI, following the nomenclature of ISteidel et al.l (120001 ) and I Weijmans et al. 

(hoioh . 



In this paper we present new spectroscopic [O ill] observations of LABI. We focus on 
the two Lya subregions C15 and Cll, the two parts of the larger LABI structure in which 
we detected [Oin]. The use of [Om] to study the kinematics of LABI is powerful because 
[O ill] is not subject to resonant scattering as Lya is. Comparing [O ill] to Lya emission 
allows us to characterize any systemic offsets between Lya and [O ill] , and thereby detect 
the presence of ou tflows or inflows. We fi rst demonstrated the efficacy of this method in a 
sample of L AEs in iMcLinden et al.l (120111 ) . 



Our new [O Hi] d a ta presented her e is co mpared to Lya data originally presented in 
Weijmans et al.l (120101 ) . I Weijmans et al.l (120101 ) measured Lya from each of the 5 subregions 
in LABI with the integral field spectrogr aph SAURON over 23.5 hours (including 9 hours 
of SAURON data originally obtained by iBower et al.l (120041 ) ). We obtained the reduced 
datacube produced from these observations from Weijmans. This allows us to extract Lya 
profiles at locations corresponding to our NIR [ Qui] observation s for c areful comparison of 
[Om] and Lya redshifts. Henceforth we refer to lWeijmans et al.l (120101 ) as W10. 



In Sections [2] and [3] we present our [O ill] detections from two near-infrared (NIR) 
spectrographs (NIRSPEC and LUCIFER). In Section [5] we look for any velocity offsets 
between our measured [O ill] redshifts and those of Lya to look for any evidence of inflows 
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or outflows and in Section [6] we discuss the implications of our findings of velocity offsets 
(Av) that are modestly negative or consistent with zero. We also compare our results to 
other authors and explore if there are currently any Lya radiative transfer models that can 
match our results. Where relevant, we ad opt the standard co smological parameters H = 70 
km s _1 Mpc -1 , Q m = 0.3, and = 0.7 ( jSpergel et al.l 120071 ). We use the following vacuum 
wavelengths, 1215.67 A for Lya, 3727.092/3729.875 A for [On], 4862.683 A for H/3, and 
4960.295/5008.240 for [O ill] from the Atomic Line List v2.O40. 



2. NIRSPEC DATA AND REDUCTION 

We ini tially detected JO ill ] emission from LABI using the near- infrared spectrograph 



NIRSPEC (IMcLean et alJll998h on the 10m Keck II telescope on 6 August 2010 (UT). We 
used the low-resolution mode of NIRSPEC, with the 42 x 0.76" slit and the NIRSPEC-6 
filter which covers 1.558 - 2.315 /im. The spectral resolution of NIRSPEC in this setup is 
~ 290 km s _1 . This filter encompasses the redshifted (z ~ 3.1) [O ill] doublet and covers 
redshifted H/3 as well, though we did not see this line. We completed three 500-second 
integrations, for a tot al exposure of 25 m inutes. The longslit was oriented so that LABI 



regions C15 and Cll ( ISteidel et al.l 120000 both lie directly on the slit. Region R2 also has 
some peripheral coverage, though not directly through the location of its peak Lya emission. 
Due to the short length of slit we were unable to place an additional continuum-bright object 
on the slit, so LABI was acquired via blind-acquisition from a nearby star. 

[O ill] detections from C15 and Cll are evident in single, raw 500-second exposures, 
when a second frame is subtracted from the frame of interest. We find these detections at 
their expected locations in the spatial direction along the slit, and find they also have the 
separation from one another that we expect for Cll and C15. In addition, the detections are 
in the wavelength range expected of each region's Lya redshift. The dither pattern we used 
is also clearly visible in positive and negative detections when we perform this sort of quick- 
look subtraction, providing assurance that these detections are not transient cosmic rays in 
a single exposure. These facts combined give us confidence that the emission we detect is, 
in fact, from [O ill] emission from C15 and Cll. See the top panel of Figure [2] for these 
detections in a sky-subtracted 2D frame. The [O ill] emission from C15 is strong and comes 
through as such in both our 2D and ID reduction processes detailed below. The [O ill] 
detections from Cll appear much fainter and are not as evident, though still marginally 
detected, throughout our 2D and ID reduction processes. 



http://www.pa.uky.edu/ ^peter/atomic/index.html 
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Initial data reduction of each 500-second exposure was done using NIRSPEC_REDUCE, 
a set of IDL programs written by G. D. Becker specifically for reduction of NIRSPEC longslit 
data. We used these scripts for flat fielding and sky subtraction in each exposure. The 
sky subtracti on process in NIRSPECLREDUCE is based on the algorithm of D. Kelson 



(lKelsonll2003l ) which provides excellent sky subtraction of even tilted skylines, such as those 
in NIRSPEC longslit data. The outcome of these reduction steps are three individual, 
2D, flat-fielded, sky-subtracted exposures. We also reduced an argon lamp exposure and a 
standard star exposure in this same way. 

We then fed these exposures into IRAF procedures in the WMKONSPEC packag^E 
Each frame was rectified to a horizontal-vertical grid in x- y using the tasks XDISTCOR 
AND YDISTCOR, which remove x and y distortion in the images, respectively. Once the 
exposures were rectified, we used IMCOMBINE to median combine the frames into a single 
exposure. We specified offsets in the IMCOMBINE procedure to remove dithers along the 
slit that were performed during our observations. The result, what we call our reduced 2D- 
spectrum, is shown in Figure [2j The locations of our [Om] detections in Cll and C15 are 
circled. 

To extract ID spectra we used the IRAF DOSLIT procedure. We first defined an 
aperture trace using a bright standard star observation which can be easily traced along the 
entire slit. We then transferred this aperture to the correct spatial locations in our science 
exposure to extract spectra of C15 and Cll. We used this transferred aperture since neither 
region has continuum emission that we are able to trace for aperture creation. The spectra 
were wavelength calibrated using an argon lamp. The RMS error from wavelength calibration 
was 0.67 and 0.45 A for C15 and Cll respectively. We used the IRAF WMKONSPEC task 
SKYINTERP to remove remaining residuals from sky lines. The resulting ID spectra from 
C15 and Cll are shown in Figure [21 

We flux calibrate our NIRSPEC spectra using a magnitude 8.85 (V band) A1V star 
that was observed in the same setup as our science observations (0.76" slit). We averaged an 
AOV and an A2V Pickles model spectrum to a pproximate an A1V spectrum. We then scaled 



down the A1V stellar spectrum (jPicklea Il998l ) to match the magnitude of the observed star. 



We created a sensitivity function with units of erg cm -2 A^ 1 counts -1 by dividing the model 
spectrum by the observed stellar spectrum and multiplying by the length of the observation. 
This sensitivity function is then multiplied by the extracted NIRSPEC spectra for C15 and 
Cll, and then the result is divided by the integration time for each object to produce flux 
calibrated spectra in units of erg s -1 cm~2 A -1 . 



2 http://www2. keck, hawaii.edu/inst/nirspec/wmkonspec. html 
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Fig. 1.— Orientation of longslit for NIRSPEC and LUCIFER observations. NIRSPEC slit 
is in white, LUCIFER slit is magenta. Boxes indicate locations of subregions C11,C15, and 
R1-R3, from W10. Im age from STIS-MIRVIS observations, retrieved from HST archives 
(jChapman et al.ll2003l ). 
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Angstroms Angstroms 

Fig. 2. — Top image is median combined, sky-subtracted, distortion corrected 2D spectrum 
from NIRSPEC. Top image has been smoothed with a Gaussian function with radius of 
three pixels. Emission lines from Cll and C15 are circled. [Om] doublet from C15 are 
upper circles, 5008.24 A line from Cll is lower circle. Bottom row contains ID NIRSPEC 
spectra for C15 and Cll (top panel) and extracted background. Best-fit Gaussians are 
overlaid on emission lines - see Section I5TT1 Blue is 4960.295 A line and red is 5008.24 A line. 
Extracted background has been multiplied by -1 so one can easily recognize the sky lines, 
but the sky lines were actually negative in some places due to imperfect sky subtraction in 
the 2D image. 
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LUCIFER DATA AND REDUCTION 



We subsequently made additional NIR observations of LABI using LUCIFER (LBT 
NIR Spectrograph Utility with Camera and Integral-Field U nit for Extragalactic Research) 
on the 8.4m LBT (ISeifert et al.ll2003l ; lAgeorges et al.ll2010l ). We used the longslit mode of 
LUCIFER with a 1" slit utilizing the H+K grating with 200 lines/mm and the N1.8 camera. 
LUCIFER has a spectral resolution of ~ 318 km s _1 near 1.6 /im and ~ 233 km s _1 near 
2.2 /im. We combine three 300-second integrations, for a total exposure of 15 minutes. We 
placed the longslit at a slightly different orientation than our NIRSPEC observations in 
hopes of capturing more emission from R2 and C15. For the LUCIFER setup, C15 and R2 
lay directly on the slit, with part of Cll on the edge of the slit (see Figure [I]). Given the 
length of the LUCIFER longslit, we were able to place a galaxy with continuum on the slit 




Fig. 3. — Third NIRSPEC frame subtracted from second NIRSPEC frame, before x-axis and 
y-axis distortion correction, and sky subtraction. [O ill] emission (5008.24 A line) from Cll 
and C15 is more evident than in Figure [2j Emission from C15 in cyan circles, emission from 
Cll in green circles. Emission from second frame is black (positive), emission from third 
frame is white (negative). Positive-negative dither pattern is clear, showing a detection in 
both frames displayed here. 
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as well. This aids in aperture extraction during the reduction process. 

We reduced the LUCIFER data in a very similar manner as the NIRSPEC data, but we 
used a modified version of the NIRSPEC_REDUCE package to accommodate the different 
detector size and orientation of the LUCIFER data. After the NIRSPEC_REDUCE proce- 
dures, the individual exposures were again median combined with IMCOMBINE and offsets 
from dithering along the slit were accounted for. Figure H] shows the combined 2D spectrum 
after this step. For the ID extraction, we created an aperture trace using the continuum 
source that shared the slit with our science targets, instead of the standard star as in our 
NIRSPEC procedure. Then we shifted the aperture to the correct spatial locations to ex- 
tract ID spectra for C15, Cll and R2. The spectra were again wavelength calibrated with 
an argon lamp exposure. Figure H] shows the ID extraction of C15 from our LUCIFER data, 
which was the only region from which we detected [O ill] in our LUCIFER observations. We 
flux calibrate our LUCIFER spectra using a magnitude 6.16 (V band) A5V star that was 
observed in the same setup as our science observations (1" slit). We used a Pickles A5V 
model spectrum and otherwise the calibration process is otherwise the same as described in 
Section [2J 



C15- LUCIFER 




2.02x10* 2.03x10* 2.04x10* 2.05x10" 2.06X10 4 2.07x10* 

Angstroms 



Fig. 4. — Left plot is 2D LUCIFER spectrum of C15. The image has been smoothed by 
a Gaussian with a radius of 3 pixels and the [Om] detection (5008.24 A) is circled. Right 
top panel is extracted ID LUCIFER spectrum of C15. Feature at ~ 20355 A is a bad 
column. Best-fit Gaussian is overlaid on emission line - see Section I5.ll Bottom right panel 
is extracted background - only minimal sky line residuals remain because sky subtraction in 
2D spectrum was excellent. Extracted background had been multiplied by -1. 
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4. Lya DATA 

We extracted Lya profiles from the reduced SAURON datacube, data that was initially 
presented in W10. The SAURON data was registered to the R band image that was used to 
align our NIRSPEC and LUCIFER slits. Next, we laid our LUCIFER and NIRSPEC slits 
(Figure [Q) on the registered SAURON data and extracted the spaxels in the SAURON data 
that corresponded with the LUCIFER and NIRSPEC [O ill] detections. This is essentially 
creating a virtual longslit aperture for the SAURON data. The LUCIFER longslit was 1" 
wide and the NIRSPEC longslit was 0.76" wide, so our virtual longslit for the SAURON 
data was 3 spaxels across, or 1.2" to provide an approximate match to the LUCIFER and 
NIRSPEC apertures. The extracted spaxels were summed to create ID spectra for Cll and 
C15. Because we observed C15 with both NIRSPEC and LUCIFER, with slightly different 
slit orientations, we extracted two spectra for C15 - one that corresponds to our NIRSPEC 
slit orientation, and one that corresponds to our LUCIFER slit orientation. When discussing 
results for C15 below we specify whether the results are for C15 from the NIRSPEC alignment 
or C15 from the LUCIFER alignment. Cll was only detected with NIRSPEC, so we only 
extracted one spectrum for Cll. The Lya lines detected in the ID spectra described above 
are shown in Figure [5j 

To get error bars for the flux in the extracted ID spectra we selected 1000 random spaxels 
in the SAURON datacube that were outside the area of the blob. At each wavelength we 
created a histogram of the 1000 flux values found in our random spaxels at that wavelength, 
fit a Gaussian to that histogram, and adopted the sigma of that Gaussian as the 1 sigma 
error on the flux at that wavelength. 



RESULTS 



5.1. [Om] Redshifts 



As we did in iMcLinden et al.l ( 120111 ) we fit detected [O ill] lines with a single symmetric 
Gaussian plus a constant, using the IDL routine MPFITEXPE0. We fit the NIRSPEC and 
LUCIFER spectra independently. The central wavelength of the best fit Gaussian deter- 
mines the [Om] redshift. We fit the 4960.295 A and 5008.24 A lines independently for C15 
in NIRSPEC and we find only the 5008.24 A line in the LUCIFER data for this galaxy. 
Given the redshift of this region from the 5008.24 A line in LUCIFER, the 4960.295 A line 



3 http: / /cow. physics. wisc.edu/ craigm/idl/ down/ mpfitexpr.pro 
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should fall at ~ 20332 .3 A, right on the edge of the 20339.497 A (vacuum) OH emission line 



(IRousselot et al.lll999[ ). This may explain why, after sky interpolation we are unable to detect 
this line in the slightly lower resolution of LUCIFER. The agreement between the redshift 
derived from 4960.295 and 5008.24 A lines in NIRSPEC spectrum is good (see Tabled]). We 
take the average of the 4960.295 A and 5008.24 A redshift as the derived systemic redshift 
for C15 from NIRSPEC, and use the redshift of the single line for C15 from LUCIFER. 

Only the stronger 5008.24 A line was detected in Cll, and only in the NIRSPEC 
spectrum, so the 5008.24 A line alone defined the redshift for this region. As mentioned in 
Section [3] the location of longslit in the LUCIFER setup was optimized for detection of R2 
and C15, so it is not surprising that we did not have a detection for Cll in the LUCIFER 
data. All these [O ill] redshifts were corrected for the earth's motion using topocentric radial 
velocities^ appropriate for the date and location of the observations. The error bars on the 
redshift are a compilation of the RMS from wavelength calibration during data reduction, 
the 1 sigma error on the best-fit central wavelength from Gaussian fitting, averaging of two 
redshifts when applicable, and a 0.02 km s _1 uncertainty on the topocentric radial velocities. 
See Table [T] for a summary of this data. 



5.2. Lya Redshifts 

To determine the redshi ft of the Lyq emission line, we followed the same method 



ology we previously used in iMcLinden et al.l (120111 ). Namely, the Lya profiles for Cll 



and C15 were fit with a single asymmetric Gaussian plus constant, using IDL routine 
ARM_ASYMGAUSSFIT§. These fits are overlaid on the extracted ID Lya spectra in Figure 
[5J The asymmetric Gaussian fitting routine allows for, but does not require, that the fit be 
asymmetric. The central wavelength of the best-fit asymmetric Gaussian defines the Lya 
redshift. 

The error on the line center was determined from 1000 Monte Carlo iterations. In each 
iteration the flux at each wavelength was altered by a random amount proportional to the 
error bar at that point and the altered spectrum was re-fit with an asymmetric Gaussian. 
The standard deviation of the 1000 best-fit line centers is used as the la error on the Lja 
line center. 

We note that the SAURON IFU datacube comes from a combination of data obtained 



4 http://fuse. pha.jhu.edu/support/tools/vlsr. html 
5 http://hubble.as.arizona.edu/idl/arm/arm^asymgaussfit.pro 
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in 2002 (IBower et al.l 12004 ) and 2006 (jWeijmans et al.ll2010l ). The data were not corrected 
to the local standard of rest (LSR) before combination into a single datacube. Therefore, 
the Lya profiles derived from the combined data may be slightly broadened with a slightly 
shifted center compared to the intrinsic Lya profile. To account for this effect, we calculated 
the corrections that should be applied to the 2002 (~ -26.6 km s -1 ) and the 2006 data 
(~ +2.1 km s -1 )B Combining two such Gaussians, offset from one another by ~ 28.7 km 
s , means that the intrinsic line center (if they had both been calibrated to the LSR before 
combination) is shifted by ~ 12.3 km s -1 . We therefore fold this 12.3 km s _1 term into 
our error calculations for redshift of the Lya line and the subsequently derived [O lll]-Lya 
velocity offsets. Table [1] summarizes our Lya and [O ill] redshift information. 



C15 (matched to LUCIFER) 




C15 (matched to NIRSPEC) 



4950 4960 4970 4980 4990 5000 5010 
Angstroms 
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Angstroms 
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Fig. 5. — Best fit Gaussians for CI 5 (top row) and Cll (bottom row) from 
ARM_ASYMGAUSSFIT show in blue. Lya spectra are shown in black. Bottom right 
panel shows Cll simultaneously fit with a double Gaussian using MPFITEXPR - the two 
Gaussians are shown in blue and red, and their sum is shown in orange. 



http://fuse.pha.jhu.edu/support/tools/vlsr.html 
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5.3. Velocity Offsets 

We derive a velocity offset ( Av) between the [0 ill] and Lya lines by comparing the 
redshifts for each line as derived from the central wavelength of the best-fit Gaussians to the 
two lines, where the fitting procedure is described above. We consider the [O ill] redshift 
to be the systemic redshift. This method finds an offset between [Om] and Lyct in C15 of 
-72.4 ± 41.6 km s~ l from the LUCIFER data and -51.3 ± 42.1 km s" 1 from the NIRSPEC 
data. These results are remarkably consistent with one another, in spite of being measured 
in slightly different apertures with different instruments. These results are modestly negative 
though also very nearly consistent with zero. We find in Av is 5.8 ± 32.9 km s _1 for Cll 
in the NIRSPEC data, which is consistent with no offset. See Tabled] for a compilation of 
these results. 

One may argue that the rather broad Lya line in Cll may be better fit with a double 
Gaussian profile, especially if one considers that the bump to the left of the highest peak is 
not noise, but in fact evidence of a second, unresolved peak. This could be a blue bump that 
is not fully resolved and sepa rated from the main re d peak. Or it could be the smaller red 



peak at v=0, as in Fig. 12 in lVerhamme et al.l (120061 ) . To consider these possibilities, we fit 



Cll a second time, simultaneously fitting two Gaussians plus a constant (see bottom right 
panel of Figure [5]). When we do this we find that the right peak yields only a modest offset, 
where the Lya line is offset by ~ 18 km s _1 from the [Om] redshift of Cll from NIRSPEC. 
The bluer Lya peak is blueshifted with respect to [Om] by ~ 808 km s _1 . This leads to a 
rather inexplicably large offset between the two Lya peaks, ~ 826 km s^ 1 , especially when 
the red peak is so mildly offset, which may disfavor this secondary interpretation of the 
modest bump as a blue bump. While acknowledgement of the additional fit described above 
is worthwhile, the discussion stays much the same, the observed velocity offset between [O ill] 
(systemic) and Lya is at best modestly negative and/or may be consistent with zero. 



5.4. [O ill] Flux in C15 and R2 

As previously described in Section I5TT} the [O ill] lines are fit with a symmetric Gaussian 
plus a constant. The resulting area under the Gaussian gives us a line flux measurement for 
each line. We find an [Om] flux of 1.1 ± 0.1 x 10^ 16 erg s" 1 cm" 2 for C15 from LUCIFER 
(5008 A line only) . For C15 in NIRSPEC we find line fluxes of 1.6 ± 0.5 x 10~ 17 erg s" 1 
cm" 2 (4960 A) and 4.6 ± 0.5 x 10~ 17 erg s" 1 cm- 2 (5008 A). Finally, for Cll from NIRSPEC 
we measure a line flux of 1.6 ± 0.4 x 10~ 17 (5008 A line only). 

Since we do not detect an [O ill] line in R2 we instead measure an [O ill] line flux upper 
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limit of < 3.4 x 10~ 17 erg s _1 cm -2 . LUCIFER has sufficient wavelength coverage to place 
upper limits on both [0 ii] and H/3 in our C15 spectrum. We derive upper limits for [0 n] 
and H/3 in C15 of 1.3 x 1CT 17 and 2.2 x 1CT 16 erg s _1 cm -2 from our LUCIFER spectrum. 
This combination of a measured [Olll] line flux and upper limits on [On], and H/3 line flux 
in C15 allow us to put constraints on the metallicity of the LBG embedded in C15 using 

R 23 , where R 23 is defined as ^° "^{P — dPaeel et all ll~979l : iKewlev fc Dopital l2002h . For 
C15, R23 > 0.76 given the line flux measurements and limits described above. NIRSPEC 
wavelength coverage only allows us to place upper limits on H/3. We find IL3 upper limits 
of 3.7 and 7.2 x 10~ 18 erg s _1 cm -2 for C15 and Cll, respectively. 

To compute the 3a line flux upper limit for R2 quoted above, we added a mock Gaussian 



emissio n line to the spectra to represent [O ill] , similar to the procedure in iFinkelstein et al. 



( 2011bl ). The sigm a of the Gaussian was fixed to 5.52 A, or the a from our faintest [Olll] 
detection to date ( iMcLinden et al.l 1201 ll ). The area under the mock line was measured 
using a symmetric Gaussian, this area determines the line flux of the mock line. Then we 
determined the noise on the line flux measurement from 1000 Monte Carlo iterations, where 
the flux array was modified each time by a random amount proportional to the error bars. 
We repeated this process, each time decreasing the area under the mock Gaussian until the 
signal to noise ratio (SNR) dropped below 5a. The line flux in the mock line where the SNR 
crossed below 5a determines a. However, because one cannot know, without an nebular 
emission line measurement for reference, exactly how much, if any the Lya line is offset from 
the [O ill] line, we repeat this calculation, fixing the mock line at different redshifts to mimic 
different velocity offsets. We found the 3a [O ill] line flux detection limit as an average of 
this technique from 6 different redshifts corresponding to velocity offsets of 0-500 km s _1 , in 
increments of 100 km s _1 . A range of - 500 km s _1 was chosen to mirror the magnitude 
of Lya - [Olll] velocity offsets we have observed of 52 - 342 km s" 1 in three z ~ 3.1 LAEs. 



Table 1. Comparison of [O ill] and Lya redshifts 



Region 


4960.295 


5008.24 


[Om] 

"Z-avg 


Zcorr 


Lyo 

z 


Av (km s x ) 


C15 (LUCIFER) 
C15 (NIRSPEC) 
Cll (NIRSPEC) 


3.0973 


3.0983 
3.0980 
3.1002 


3.0976 


3.0986 ± 0.0003494 
3.0974 ± 0.0005095 
3.0999 ± 0.0003590 


3.0976 ± 0.0004159 
3.0967 ± 0.0002070 
3.1000 ± 0.0002138 


-72.4 ± 41.6 
-51.3 ± 42.1 
5.8 ± 32.9 



Note. — Comparison of [Om] and Ly« redshifts in Cll and C15 as measured with LUCIFER and NIRSPEC. z aV g is listed 
when the redshifts from 4960.295 and 5008.24 A lines were averaged. z corr is the redshift after correction to the LSR. Av is 
velocity offset between [Om] and Lyes. 
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The range of 3a detection limits over this wavelength range was 3.12 - 3.54 x 10 -17 erg s -1 
cm -2 . The [Oil] and H/3 upper limits described above are found using the same procedure, 
except the upper limits need only be derived at a single fixed wavelength for each line, the 
wavelength defined by the [Om] redshift. 

Given [O ill] line flux detections in C15 and Cll and an [O ill] upper limit in R2, we can 
compare the nebular emission from these subregions. This is of interest because C15 and 
Cll are associated with LBGs embedded within the larger LABI halo structure, whereas 
R2, in spite of its str onger Lya emission , is not associated with any u nderlying galaxy 



fWe ijmans et al.ll2010f ). The ratio of Lya luminosity in C15 to R2 is 0.7 ( IWeijmans et al. 
2010|). We find that the ratio of [Om] in C15 (LUCIFER) to R2 is > 3.4 and the ratio of 
[Om] in C15 (NIRSPEC) to R2 is > 1.4, meaning that while R2 is brighter in Lyct C15 is 
brighter when looking at [O ill] nebular emission. Finally, the [O ill] to Lyct ratio measured 
in the region without an LBG is therefore at least 1.9 - 4.5 times smaller than the same 
ratio measured in the LBG. This measurement would likely indicate that something other 
than star formation is powering the Lya emission in region R2 and that there may very well 
be different sources powering Lya emission in different regions of the same blob. 

Cll presents a somew hat different story. The ratio of Lya luminosity in Cll to R2 is 



1.0 flWeiimans et alJl2010f >. Then the ratio of [Om] in Cll (from NIRSPEC) to R2 is only 
> 0.5. The [Om] to Lyct in Cll is only 0.06, which is smaller than that same measurement 
in R2, where the ratio is < 0.13. These divergent results are interesting, perhaps giving 
further credence to the idea that there are different processes at work in different regions of 
this blob. 



5.5. Asymmetry and [O ill] - Lya Offset 
We have now mea sured the velocity offset between [Om] and Lya in 3 z ~ 3.1 LAEs 



(IMcLinden et al.l 120111 ) and added 3 Av measurements for LABI in this work. Amongst 
these measurements, the regions Cll and C15 in LABI have the smallest absolute Av mea- 
surements and are the only negative values. This is in contrast to three z ~ 3.1 LAEs in 
which we found offsets ranging from 52 - 342 km s _1 . Given this information, we can com- 
pare another signature of outflows, namely asymmetry in the Lya profile, with the velocity 
offset measurements. 

We quantify asymmetry as a r ed/&biue where a re d is the sigma of the red side of the 
best-fit asymmetric Gaussian and ouue is the sigma on the blue side of the best-fit asymmet- 
ric Gaussian, where both are parameters returned by the ARM_ASYMGAUSSFIT routine. 
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C15 ([OMI] from LUCIFER) C15 ([ONI] from NIRSPEC) 




C11 ([ONI] from NIRSPEC) 




Fig. 6. — Plot of Lya (blue) over [O ill] (black) where Lya has been shifted to [O ill] frame 
via Lya:rQ in j = Lya k serve( j x ffjff^- Observed spectra are shown with solid lines, best 
fit Gaussians for both Lya and [Olll] are shown with dashed lines. Top panel is C15, left 
plot shows [O ill] from LUCIFER, right plot shows [O ill] from NIRSPEC. Bottom panel is 
Cll from NIRSPEC. 
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With this definition, a profile with asymmetry > 1.0 is considered asymmetric, and the 
asymmetry is dominated by the red-wing. Objects with asymmetry = 1.0 are symmet- 
ric, < 1.0 have blue- wing dominated asymmetry. The red-wing dominated asymmetry is 
the expected direction of the asymmetry in the Lya line from high-z gal axies, as the red 



side o f the line can be enhanc ed in the presence of an expanding shell (jVerhamme et al. 



2006 



( iRhoads et al 



2003 



Diikstra fc Wvithei 120101) and/or b y interaction with neutral Hydrogen in the IGM 



Dawson et all 120041 ) . Measured in this way Cll has an asymmetry 



of 0.74 ± 0.02 and C15 has an asymmetry of 0.9 ± 0.03 and 1.1 ± 0.03 (LUCIFER and 
NIRSPEC apertures, respectively). An asymmetry of < 1.0, like Cll, can suggest infalling 
material. However, the Lya-[Olll] velocity offset of 5.8 ± 32.9 km/s in Cll is not very 
supportive of such an interpretation. In three z ~ 3.1 LAEs in which we have measured a 
velocity offset between [Olll] and Lya we find asymmetries of 0.97 ± 0.1, 1.04 ± 0.1 , and 
1.65 ± 0.1 for LAE7745, LAE27878 and LAE40844, respectively iMcLinden etaPboilh. W e 
add fo u r additional asymmetry data points by using four z ~ 2 LAEs from lHashimoto et al. 



(120121 ). lHashimoto et all (120121 ) measured redshifted Lya lines in these objects with respect 
to Ha lines in the same objects. They report velocity offsets of 18 - 190 km s _1 in their 
four LAEs, in good agreement with our range of 52 - 342 km s _1 in LAEs. We measured 



the asymmetry of the Lya lines presented by lHashimoto et all (120121 ) in the same man- 
ner as above, by fitting each line with asymmetr ic Gaussian and quantif ying asymmetry as 
Cred/cbiue- Measured in this way the LAEs from lHashimoto et all (120121 ) have asymmetries 
of 0.53 - 1.7. 

Figure [7] demonstrates a possible trend where asymmetry in the Lya profile increases 
with increasing velocity offset. We find a moderate Pearson linear correlation coefficient 
of 0.427 (P=0.110) between the velocity and asymmetry values, suggesting the trend may 
be real. This velocity-asymmetry correlation is not unexpected, as increasing (red-wing 
dominated) asymmetry in the Lya profile is tied to increasing shell expansion velocities in 
Verhamme et all (120061 ). when outflows are modeled with a central monochromatic source 
and an expanding shell. This is because Lya photons are seen as more redshifted by the 
neutral hydrog en in the shell with high er shell expansion speeds, decreasing the cross-section 
for interaction (jVerhamme et al.ll2006l ). A symmetric line, as seen in C15, is not expected in 
models with large outflows, but could be consistent with static/near ly-static profiles - if the 
two symmetric peaks produced from a static slab or shell scenario (jVerhamme et al.ll2006l ) 
are unresolved. The authors note that the point attributed to LAE40844 strongly influences 
this correlation. Removal of this point and recalculation of the Pearson coefficient with the 
remaining eight points only yields a coefficient of 0.170 (P=0.331), which does not suggest 
any significant correlation. 
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DISCUSSION 



6.1. Comparison with Blobs at z ~ 2 



We have measured A v between [0 ill] and Lya in two subregions of LAB 1 and found 
that Av is modestly negative or consistent with km s" 1 . These measurements seem to 
downplay the role of powerful outflows in powering the Lya in these regions and could even 
hint at some infall. Interestingly, this measurement of Av = small is not the first time this 
phenomenon has been reported in a Lya blob, suggesting this is an impo rtant phenomenon 
that must be explored to better understand the nature of high-z LABs. I Yang et al.l (120111 ) 
investigated two z ~ 2.3 LABs, where one blob (called CDFS-LAB02) had Av ~ 230 km s" 1 
and the oth er blob had Av co nsistent with zero (called CDFS-LAB01). The velocity offsets 
reported in lYang et al.l ( 1201 ll ) were measured by comparing the redshift of Ho to that of 
Lya, a very similar tactic to our comparison of [Om] to Lya, as both [Om] and Ha are 
nebular emission lines and therefore probe the same regions. 



6.2. Comparison to LAEs and LBGs 



Our Av measurements, and those of I Yang et al.l ( 120111 ) , are significantly less than the 
larger velocity offsets typically seen in Lyman break galaxies (LBGs) and even Lya emitting 
galaxies (LAEs), velocit y offsets which are t ypically interpreted as clear signatures of strong 
winds in these galaxies. ISteidel et al.l (120101 ) report median velocity offsets between Ha and 
Lya of 445 km s _1 in 41 z ~ 2.3 LBGs. And even LAEs, whose typical velocity offsets have 
been found to be smaller, have Av as large as 342 km s~ 1 (jMcLinden et al.ll201ll ). 



Steidel et al.l (120101 ) also report a median velocity offset between Ha and strong inter- 
stellar absorption lines of -164 km s -1 in 86 LBGs, which further supports interpretations 
of the presence of outflows in LBGs, si nce the blue-sh i fted a bsorption implies absorption in 
material moving towards the observer. IShapley et al.l (120061 ) previously measured a redshift 
from low-ionization interstellar absorption (LIS) lines in the LBG Cll, finding that the ab- 
sorpt ion lines are offset fr om the Lya line by -380 km s _1 . Using the LIS redshift of 3.0962 
from IShapley et al.l (120061 ) and comparing this to the Lya redshift we derive for Cll yields 
an offset of ~ 278 km s _1 . Comparing our [Om] (systemic) redshift for Cll (3.0999) to 
the LIS redshift, we can estimate that the LIS lines in Cll are offset from [Om] by ~ -270 
km s -1 . This comparison of Lya and [Om] redshifts to LIS absorption redshifts provides a 
stronger signature of an outflow than we get when comparing Lya and [Om]. In fact, it is 
particularly interesting that the magnitude of this second signature of winds (i.e. blueshifted 
interstellar absorption lines) is so similar when comparing the Steidel LBGs to Cll, yet the 
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magnitude of the offset between nebular emission lines (Ha or [Olll]) and Lyct offset are so 
different. We note that the lack of Av between [0 ill] and Lyct does not have to rule out 
some outflows in LABI. Rather, the lack of a Lya-[0 ill] offset may just imply that outflows 
are not a significant mechanism for helping Lyct photons escape. 



6.3. Previous Studies of LABI 



This phenomena, i.e. Av = small and/or = km s" 1 , leads to the question of whether 
the lack of substantial velocity offset between Ha-Lya or [0 m]-Lya in these blobs in fact 
rules out outflows or if there is some, yet to be understood phenomena, that damps or erases 
this particular wind sig nature. This question is particularly relevant given the recent report 
from lHayes et al.l ( 120111 ) that there is polarized radiation emanating from LABI, polarization 
that is indicative of scattering of Lya photons at large radii with respect to their site of 
production. This may be a sign of outflows helping to drive the Lya photons to these large 
radii, but we are not seeing the velocity offsets between nebular emission lines and Lya that 
we would expect to see if this was the case, velocity offsets we have been able to see in other 
obj ects at similar redshifts believed to have strong winds. In addition, iBower et al.l (120041 ) 
and IWeij mans et al.l (120101 ) both measure a velocity shear in the Lya emission from Cll and 
C15. As the authors point out, such a velocity shear could be consistent with infalling gas, 
outflowing gas and/or rotation of the system, and such scenarios cannot be differentiated 
from the Lya data alone. While both papers use this velocity shear to argue in favor of 
the presence of outflows in Cll and C15, we can report no signature of such outflows when 
we compare the redshifts of [Olll] and Lya, a comparison that has proven to be a useful 
diagnostic of winds in LAEs and LBGs at similar redshifts. 



6.4. Comparison to Radiative Transfer Models 



We explored available Lya radiative transfer models to see if there were any models that 
might shed light on the physical conditions that could lead to a modestly negative Av or Av 
of ~ km s" 1 be t ween [O ill] (or Ha) and Ly<x We focused on the Lya profiles produced by 
Verhamme et al.l ( 120061 ) (henceforth V06) from their 3D Monte Carlo Lya radiative transfer 
code. V06 explore a variety of physical conditions and geometric orientations to explore 
the diversity of Lyct profiles that arise from different conditions. See V06 for extensive 
details on these models and the model parameters. We find, however, that none of the 
m odels presented in V06 are in great agreement with our observations (or the observations 
of lYang et al.l ( 1201 ll )). The only models that are marginally consistent with our observations 
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are those that have two significant Lya peaks, where the centroid of those two peaks is 
centered at v ~ 0, with one peak redward of v= 0, and the other blueward. Convolved 
with an appropriate Gaussian to mimic the instrumental resolution of SAURON (290 km 
s _1 ), these models somewhat resemble the Lya profiles presented here. This double-peaked, 
centered at v ~ profile occurs when a central monochromatic source is embedded in a 
static slab (with or without dust, V06 Figure 3), a central monochromatic source sits in a 
non-expanding shell (V06 Figure 14), and when a single monochromatic source sits in an 
expanding shell with very small velocity gradient (Vmax = 20 km s _1 , V06 Figure 7). The 
closest match (see Figure [8]) is from V06, Figure 14 with a static shell. In addition, we 
explored whether the infalling halo models of V06 might agree with our observations, these 
models are presented in Figure [91 There may be some qualitative agreement between the 
models and observed profiles here, but the models explored are too far blueward of v=0 to 
match the Lya profiles we present here. In either scenario (expanding or infalling) the main 
peak is expected to be asymmetric, which is inconsistent with the Lya profiles of at least 
C15 presented here. We contend that more modeling of Lya radiative transfer with direct 
application s to th e observations we have presented here and those observations presented in 



Yang et al.l (120111 ) needs to be done, to better understand the physical conditions, geometry, 



and kinematics that can produce single peaked Lya lines, with Av values that are modestly 
negative or consistent with zero. 



7. CONCLUSION 

We have measured [Olll] in two subregions of LABI, Cll and C15, regions that are 
associated with underlying LBGs within the larger halo structure. We have quantified the 
velocity offset between [O ill] and Lya redshifts in these regions, finding that both mea- 
surements are modestly negative or consistent with zero. This is an intriguing result since 
powerful outflows have been proposed as possible ways to explain the luminosity and large 
spatial extent of LABs. We cannot completely rule out the presence of strong winds and 
outflows in LABI, but we can state that we do not see two typical markers of their presence. 
(1) The aforementioned result that we do not find that Lya is redshifted with respect to 
[O ill] in the two LABI subregions C15 and Cll. And (2) we do not measure strong red-wing 
dominated asymmetries in the Lya profiles of these objects, in contrast with z ~ 2 and z ~ 
3.1 LAEs where the asymmetry of the Lya line appears to increase with increasing velocity 
offset between [Olll] (or Ha) and Lya. If outflows are present in LABI, they do not appear 
to be a crucial mechanism driving Lya escape. 

In addition to the conclusion above, we have placed an upper limit on [O ill] line flux 
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from region R2, a subregion of LABI unassociated with any known galaxy and compared 
this to the [O ill] flux from subregion C15, which is associated with an LBG. We find that in 
spite of the strong Lja emission present in R2, the [O ill] flux from C15 is stronger than that 
of R2 by at least 1.4 - 2.5 times. This measurement may indicate that diverse sources of 
Lja emission may be responsible for powering different regions even within the same blob. 
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Fig. 7. — Asymmetry of the Lya profile as a function of velocity offset in five Lya-emitting 
objec ts at z ~ 3.1 (shown i n black) and four z ~ 2 objects (shown in red). Three LAEs 



from iMcLinden et al.l (120111) are labeled with the prefix 'LAE.' Red points with the prefix 
'cosmos' or 'cdfs' are from lHashimoto et al.l ( 120121 ). Overall, asymmetry may increase with 
increasing velocity offset. 
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Fig. 8. — Lya profiles (black dashed line) compared to V06 outflowing shell and halo models 
(solid blue). V06 models have been convolved with a Gaussian to match instrumental resolu- 
tion of SAURON. Top row shows comparison to model in V06 Figure 14. Bottom row shows 
comparison to model in V06 Figure 7. While these models and those in Figure [9] represent 
the closest matches to the observed profiles, clearly none are an excellent fit. 
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Fig. 9. — Lya profiles (black dashed line) compared to V06 infalling halo models (solid 
blue). V06 models have been convolved with a Gaussian to match instrumental resolution 
of SAURON. Top row shows comparison to model in left panel of Figure 5 in V06, bottom 
row shows comparison to model in right panel of Figure 5. While these models and those in 
Figure M represent the closest matches to the observed profiles, clearly none are an excellent 
fit. 



